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ABSTRACT 


An experimental investigation is presented for both steady 
and pulsatile flow of blood through the aorta. A mixture of glycerine 
and water was used as a blood model. The effect of varying the significant 
parameters, Reynolds Number, Re, unsteadiness parameter, a, and normalized 
velocity fluctuation, X, is revealed through their influence on the mass 
flow ratios; the mass flow ratios y, (n=l, Coed >) ebelag the ratios 
of the mass flows in the right subclavian, right carotid, left carotid, 
left subclavian and thoracic aorta, respectively, to the total mass flow. 
An experimental apparatus and an artificial aorta were 
designed and constructed to simulate the physiological range of these 
parameters. One type of rigid artificial aorta was examined. 
For pulsatile flow, Y, was found to be independent of the 
unsteadiness parameter, a, in the range 10 < a < 20, and was found 
to be a decreasing function in the lowest range 8 < a < 10. Furthermore, 
Y, was found to be a decreasing function of both the Reynolds Number and 
the normalized velocity fluctuation up to the range of Re = 1500. For 
larger Re, as for steady flow, Ye is a decreasing function of Re only. 
The mass flow ratio, Yos showed the same dependence as vy 
except that it was more strongly influenced by Re and not as strongly 
influenced by the normalized velocity fluctuation. The mass flow 
faci, ¥3> showed the same dependence as Yq: The mass flow ratio, 
Yas Was similarly dependent as yy except that Yq was not an inverse 
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The mass flow ratio, Yo» representing the continuation of the 
flow of blood to the body from the aortic arch was found to be 
independent of the unsteadiness parameter at least in the range 
10 < a < 20 for a given Re. Within the experimental accuracy observed 
the dependence of Y5 on X could not be evaluated. 

A visualization technique revealed the existence of four 
interdependent regions of separation for both steady and forward 
pulsatile flow. The formation of vortices was observed in the separated 
region of the main branch and double helicoidal flow was observed in 


branches four and five of the artificial aorta. 
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CHAPTER I 
INTRODUCTION 


1.1 Introductory Remarks 


An important cause of death is diseases and defects of 
the cardiovascular system. Thus fundamental research into this 
system is being pursued. The flow of blood through the vessels 
is one of the aspects of this system well suited for engineering 
research. It is essentially a problem in fluid mechanics, albeit 
a rather complicated one. 

The purpose of this thesis is to ascertain the dependance 
of the flow rate in the downstream branches of a certain artificial 
aorta on the non-dimensional parameters: Reynolds number, unsteadiness 
parameter and velocity fluctuation parameter. The high degree of 
complexity of the physical laws and equations describing blood flow 
necessitate that the study of this problem is experimental rather than 


analytical. 


1.2 The General Description of the Aorta and Its Function in 
the Blood Circulatory System 


The aorta is the great trunk artery that carries blood from 
the heart to be distributed by branch arteries through the body. The 
aorta is a wide tube with thick walls which are largely composed 
of- elastic tissue!; and like any other elastic structure, its capacity 
is determined largely by the pressure of the blood it contains. 


When blood is expelled during ventricular systole, the 
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aortic pressure rises, the aorta is distended and so accommodates a 
large part of the blood expelled, the remainder escaping through the 
arteries. During ventricular diastole, the tension in the aortic 
walls maintains the flow of blood onwards through the arteries, 
and the aorta diminishes in size until it is again distended at 
the next heart beat. In this purely passive way the aorta (and to 
a less extent its larger branches which are similar in structure) 
converts the intermittent flow from the heart into a continuous 
though pulsating flow in the arteries. As can be seen in Figure 
1.1* the blood passes from the aorta into the main branches. The 
arteries divide repeatedly until they become capillaries. (See 
Figure 1.2*). With each subdivision, the arterial lumen becomes 
smaller and the arterial wall becomes more predominantly muscular. 
By the contraction or relaxation of their muscular walls, the arte- 
ries (particularly the smaller branches or arterioles) are capable 
of great variation in calibre, and it is this calibre which largely 
determines the proportion of the total cardiac output distributed 
to the individual organs. 

Moreover, the total cardiac output could be influenced by 


e state that arterial dis- 


arterial diseases. Mitchell and Schwartz 
orders are known by such names as arteriosclerosis, atheroma, and 
atherosclerosis. The most widely adopted term, atherosclerosis, 
will be used in this thesis. It will refer to the buildup of fatty 
and/or fibrous plagues which occlude the vessel and destroy the 


a a eS 


* X-Ray pictures were taken at the University of Alberta Hospital. 


oe . oe 
§ eitiboemeas: dian sikicaeiahadiel 
vite Mguondt pniquozs ébntanin att fat i et 


i 
AL 


atevan efi at cabana’ adh <élopeatty Jatestidebeindtei 

eet setde ald Auer ebikwno boot Fo: wol¥ elt ae atfoe 

te bildaseth ntsid ef dt Thay este ah asddiqtetb aves wee: 
of bas) arioa Sis Yew’ aviesen (feud ata nt dead tyaed dixon aa! om, 

lomitaate nt rstTate ave. Hofde aauaderd “rage! ext teodae zest . a0 sie 

egauntsnos 6 cint tvesd ant jerk wor? qnedtiismnt a2 canal ’ 


suiphicnl wabeced tar 2A .zetiotze ate nt got gatesetuq Aguedt Ws 
ait .zatibvard ntom ofS ofal a2are 979 aOyt 2enkew boots oe #0 _ 
abe} .2sianliiqgs smosed Yair trem y(hemnadte. aipeth eeivete 
sainoiad adel fetve ise dt .aaretyl bay Rone deh tet ugh f 
1eluacu llr ome at ia it Tow faraedie ote ore oat tame : : 
-s$i6 ats .2tlaw seule’ sfory sg patdexetsy no Seldzersno> af? ee 
stddges Sie (eafolesdss xo esaneyd sel fame add Yirefuptoesag) wh . 
ufone! Gofiwerdti(as cin? zt Tt Shs .eidiies af notietrey seen io 0 
hotudiaserh tuytue seth yes) [ator sd) To aobinogoig ets corteretab. a 
- 2raG 4D tastiteEhnt wis ay 
yd tenet tat sc bluas sugdue oni Tetot dat rewoeelt 
nett tetra sodt aicte “tixanind the Itsitat | .aeeesrth fabystrg, 
bos. -pmcrenits (2kAneTazotarIn 26 Rsmr Move YW outed oe aNORG, 
—_ tkeret secrete unot batyaun yfubhe tet aAT .charelosernedia, a 
18? Yo qubk ted add Gt vetey Pte aT. -etasus chiy at bpew od ENRON 
eas worden tne fassan tf stihond fatiny enue sonst? sa \bep 


stetgaat atisdta. Sore mo ve nt a . 
Cylon eee 


oy 


elasticity of its walls. 


1.3 The Cardiovascular System* 


The human cardiovascular system consists of two pumps 
connected in series with two circulations as shown in Figure 1.3. 
The left ventricle pumps freshly oxygenated blood through the aorta 
into the systematic circulation which leads to arteries going to 
the head, internal organs and limbs. As the various organs are 
perfused with this oxygen-rich blood, oxygen is diffused from the 
blood into the cells of the organs. As the blood becomes depleted 
of its oxygen supply it returns to the right atrium and 
right ventricle by means of lesser and greater veins. The 
right ventricle pumps the blood through the pulmonary circulation 
(lungs) where oxygen from the atmosphere is diffused into the blood. 
The oxygenated blood then continues on to the left atrium and left 
ventricle by means of the pulmonary veins and the cycle thus continues. 
The two ventricles work in unison, both discharging their 
load of blood simultaneously into their outlets or arteries. This 
phase of the heartbeat the expulsion of the blood by contraction of 
the ventricles, is called systole. The phase in which the ventricles 
are filled with blood, each from its own atrium, is called diastole. 
The left and right ventricles must pump an equal volume of blood 


over a period of time. 


* This section comes before the literature survey in the hope that 
any unfamiliar terms used in the survey will have been explained 
adequately to those readers unfamiliar with a physiological 
terminology. 
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Figure 1.3: Scheme of the Human Cardiovascular System 
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1.4 Literature survey 


As mentioned previously the blood flow problem lies in 
a wide and complicated field. At the present time, the exact 
nature of blood flow in the human ascending aorta is not known. 
Particularly the elastic deformations of the vessel make 
a theoretical approach very difficult. Theoretical considerations 
as well as results obtained from model experiments by J.G. Porjé>, 
confirm the supposition that the effects of viscosity and elasticity 
are of secondary order of magnitude and that the equations valid 
for calculation of nonviscous flow will therefore give a reasonable 
approximation of the calculation of instantaneous flow in the ascen- 
ding aorta. Fully developed flow can only exist in the flow field 
distant from the ends of the tube where end effects are negligible. 
Atabek* considered these end effects in developing an inlet length 
theory for periodic flow in tubes with rigid walls. 

Kuchar and Ostrach” considered entrance effects associated 
with laminar flows of viscous fluids in circular elastic tubes. Of 
particular interest was their derivation of the significant dimen- 
sionless parameters. Chapter II of present. work. reviews these dimen- 
sionless parameters. 

The problem becomes even more complex in the case of 
branched vessels and, consequently, most of the work in this area 
is done with the help of experimental models rather than adopting 


a purely theoretical approach. 
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G.S. Malindzak ar? discussed the derivation of a system 
of mathematical equations to describe the transmission and reflection 
of arterial pressure pulses in the mammalian circulatory system. His 
theoretical analysis involves calculation of incident and reflected 
wave parameters derived from dependence of time, space and pressure 
data. Specifically, the paper describes the theoretical basis for 
numerical analysis resulting in values for reflected energy » phase 
velocities in the aorta, as well as experiments in which such analysis 
has been applied in the anesthetized dog. 

Robert H. Cox! reviewed the linearized wave propagation 
model which can be used to predict both transmission characteristics 
and hydraulic fluid impedance values upon system input parameters. 
Most of the models considered the propagation of small, axisymmetric, 
harmonic disturbances through an incompressible, Newtonian fluid 
contained within a distensible tube of long length. 

J.W. Amyot et a1® measured the sequential velocity develop- 
ment in the ascending aorta using a conical hot-film probe in six 
anesthetized dogs. Sequential velocity developments during systole 
are presented and instantaneous volumetric flows are obtained from 
the integrated profiles. The calculated flows are in agreement with 
standard measurements of cardiac output. The measurements show that 
maximum reverse flow occurs at the end systole when the dicrotic 
notch is observed in aortic pressure trace. Nonsymmetry was observed 
in the sequential velocity profiles during the early phases of ejection, 


followed by more nearly symmetric profiles as the flow decelerated. 


wmareye 6 Ye neni boaaunerd xt 
aeizzol wt tne nota tmenara edt sdinesh-ot -¢ pa. 18: 
ei mabey danrinayeiespratonpearenyreere 
barpeltes bw anettant Yo notistuohen: zevlount efegfans [sotseveed? 
giuazetg bane oosge -omts 10 seWatonagalt mont bovine 2ratwnwneg ovine 
aot atand [esthesoans oft’ gadirigeeh tageg ant ef fpotytoane. ab 
gently ,yerats basset te: vat aeufav at gorttoney zhayfens {sotvemun 
zraytene dove dotdw nt 2tese?nagne 26 isw 26 .atroe ant at esti toofey 
Ooh bantauritzens od? at betigge need 2ert 
nottepsqeg sysw bestrgsatt, and bows f yer \yed A Suades 
astsalhe? ‘ovate aatnabbentrt Wide soibewq of bean of hed dzivw febon 
vetoes tagnt meteez bogu! eaulen ‘sonsbignt bruit ot testi brs 
opine ive .Pitaie Yo noteapsqore att bevabtenos 2fsbue aft? to 320M 
hruft ceteetwel .nhcheesramean! ie dquewit seocedeureth ateowred 
 stgesl prof te edad eldtanetzth ¢ atittw behets - 
sddfawan’ qstroten Fetineupee ote Garuanee “fet town vi. 
xhe Wh-wdene af f¥-fod Fsxtneo '§ gotzu sd768 Ertbesoze ed? af deem 
efoteve ontwh etnaansiovet vtrofsy fatinsues2 .2gob basttemteans 
cent tpntatde orm cuolt-atiiewufoy eummatantent tas Sasnseang wre | 
dgiw Sremeniee nt eve ewott bodalwatss aaT zation batorpstat ats 


deitd wale ejremerutese onT, fugfeo sbTines Yo zseammutsen brebnste 


shirt’ ait min eloderga’ bits att te 2wor wolt seve: muatxen | 
i ceieantienthiiiaiiaiel er nt bavesedo 2b dodon 


epee fren sion. yd bewol fot 7. 
a 


ih ‘s 


Donald L. Fry et a1? derived a method of determining the 
instantaneous aortic blood velocity from the Navier-Stokes equation. 
Experimental support for the derived formulae was obtained in a flow 
model and in a living animal aorta. 


Marrill P. Spencer and A.B. Denison!? 


presented the exact 
blood flow curves of the aorta and established their validity by 
their relationship to the differential pressures existing along 

that structure. They pointed out some artifacts that may Re intro- 
duced into flow curves when recorded under certain conditions which 
considerably alter the normal physical characteristics of the arterial 
system. 

M.F. Snyder et all! devised a model of the human systematic 
arterial tree based on a lumped-parameter-circuit approximate form. 
This model has been set up and studied on an analog computer. A 
feature of this simulation is the division of the arterial system 
into sections whose lengths are inversely proportional (approximately) 
to their cross-sectional area - or what they termed "equal-volume" 
modeling. The simulated pressure and flow waveforms obtained with 
the model compare favourably with data recorded from the normal 
adult human, and exhibit such well-known features as distal delay 
and peaking of pressure pulses. The predicted aortic input impedance 
dependency on frequency checks well against measurements on the human. 
The model also provides a simple means for determination of cardiac 


output, cardiac work and cardiac power under various assumed condi- 


tions such as variation of heart rate. 
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1.5 Objectives of the Thesis 

The most sophisticated blood flow in arteries seems to be 
within the aorta. It is the purpose of this thesis to study the flow 
through the main branches of one such bifurcation and to find out 
what influences the mass flow ratios Yy> Yoo Y3> Yq and V5 throughout 
the artificial aorta, particularly when the significant dimensionless 
parameters are varied. 

In addition to the dimensionless numbers used by Rodkiewicz 


2 


and Howell ~, it was decided to incorporate the characteristic values 


introduced by 1to!3, Furthermore a modified Dean number Km (product 
of the Dean number and a constant) will be introduced. 

Some qualitative observations of flow characteristics in a 
bifurcation will be made using the visualization technique used 


by Rodkiewicz and Roussel! *, 
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CHAPTER II 
THE MATHEMATICAL MODEL 


2.1 Introductory Remarks 


The solution of the mathematical model to encompass 
and completely describe the physical processes occurring through- 
out the blood circulation system is not possible at the present 
time. The mathematical model is described in the next section, 


and is subject to the following assumptions. 


1. Blood is considered to be a homogeneous incompressible 
Newtonian fluid with a constant viscosity. 

2. The vessel will be considered to be circular. 

3. The tube walls of the bifurcation will not be flexible. 

4. The tube will be considered to be externally constrained 
in such a manner that no longitudinal displacement or motion 
of the bifurcation wall occurs. 

5. This model applies to the macrocirculation (large arteries 
and veins), but not to the microcirculation (arteries, 


capillaries, and venules). 


2.2 Derivation of Parameters 
The purpose of the mathematical models as far as this thests 
is concerned is to allow the derivation of the significant parameters. 
The problem under study may be thought to be composed of two classical 
cases of the flow, through the straight pipe and flow through curved 


pipe. The two distinct regions are shown in Figure 2.1. 
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Kuchar and Ostrach® and Kuchar and Scala? developed the 


following non-dimensional parameters for pulsating flow in region 1. 


UR 
Reynolds' Number, Re = = 


Strouhal Number, S = salsa 
UT, 
Ro 
Length Ratio, a a 
Euler Number, E = AE 
oU 
The unsteadiness parameter ,a, developed by Kuchar and Scala!? 
is as follows: 
URSteae ark 
= et vaso. = si: cum mes 
oO REexXes: xX Le ara OT, x 7 
= |? y patie (209) 
) vig 


If the pulse rate is defined as w = = equations (2.2) become 
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R =— (233) 
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Figure 2.1: The Composition of Two Distinct Regions 
of the Aortic Arch 


Figure 2.2: The Coordinate System Employed to Study 
Flow in Region 2 
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In this region, it should be noted that Re is calculated at the 
entrance to the bifurcation. 


2 


Rodkiewicz and Howell ~ considered three parameters, Re, a, and 


A where a, in equation (2.4), is a modified form of the unsteadiness Para- 
meter used by Kuchar and Seale 

Region 2, which is composed of the main arch (upper and lower 
aortic arch) of the artificial aorta, is considered to have a constant 
cross-section. Examination of several x-rays showed that though a change in 
diameter along the arch length does exist this deviation was neglected. 

The preceding considerations concerning pipe flow are valid only 
for straight pipes. In curved pipes there exists a secondary flow!6217 
because the particles near the flow axis which have higher velocity are 
acted upon by a larger centrifugal force than the slower particles near 
the walls. This leads to the emergence of a secondary flow which is directed 
outwards in the centre and inwards (i.e. towards the centre of curvature) 


near the wall. 


The characteristic dimensionless variable, which determines the 


influence of curvature in the laminar case, is the Dean Number !8219 K, 
R 
at BO! 


The coordinate system employed to study region 2 is shown in 


Figure 2.2. The z axis is perpendicular to the plane of the page. 


2.3 Analysis of Governing Equations in Region 2 


2.3.1 The Continuity Equation 


The continuity equation for cylindrical coordinates is given by 
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and introducing the following characteristic values used by H. Ito!3, 
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V5 v, : ae ee (228) 
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where, ce is average velocity, and ¢ is in range 0 < $ < 1; substi- 
tuting into the continuity equation (2.6), assuming p to be constant, 


the equation is reduced to the following non-dimensional form 
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Consider the nondimensional parameter — the leading 
coefficient of the second term of the above equation. Since s = R@ 


and R = RoR/R, this parameter can be expressed as 
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eee dae asa nou RT 
ve v v R (2.9) 
6) 
Therefore, the continuity equation can be rewritten as 
19 ~-c R 190 vo od o7o 
==(FV) +0 5-1-1 -—,) +2 (4) = 0 (2.10) 
ror, o 8 6Y OA dz 
2.3.2 The Momentum Equations 
The momentum equations governing steady flow in 
cylindrical coordinates are given by 
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Substituting (2.7) and (2.8) into momentum equations (2.11), 
and dividing the first and the second by soe. and the third equation 


by de wi. the equations are reduced to the following non-dimensional 


forms. 
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Consider the dimensionless groupings of equations (2.12). 
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which can be written as 
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Since the radial pressure gradient term cannot be neglected, choose 
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then body forces in the radial momentum equation are very smal] 


when f eee tes 
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(2217) 
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(2.19) 


then body forces in the azimuthal momentum equation are very small 


when f.. >> 1,7 17e; Le is very large and/or s is very small. 


W) 
Equations (2.16) and (2.13) show that 


Thus the final form of the non-dimensionalized momentum 


equations is 
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CHAPTER ITI 


THE EXPERIMENTAL MODEL 


The arterial system is, at its simplest, a system of con- 
duits carrying fluid from a source pump to the capillaries. That 
this system defies elegant and accurate analysis can be attributed 
largely to the fact that the pump never delivers a stroke which is 
exactly the same as any other stroke; no two aortas have exactly the 
same dimensions and mechanical properties; no system of branches ever 
exactly duplicates any other system of branches; all of the tubes are 
tapering and are made of materials which exhibit non-linear behaviour, 
and they are filled with fluid which is non-homogeneous and anomalous 


in behaviour. 


3.1 The Range of Study 


The expressions for the Reynolds number and unsteadiness 
parameter as described by equations (2.3) and (2.4) have been developed. 
It is now necessary to determine the range of these parameters for the 
physiological case. To do this, the following data has been gleaned 
from Wintrobe?, Tuttle and Schottelius<’, and Bell, Davidson, and 
Scarborough: | 

1. Blood density wauies from 1.048 to 1.066 gm/cm> 
2. The relative viscosity of blood with respect to distilled 

water varies from 3.5 to 5.4. 

3. The kinematic viscosity of distilled water equals 0.010027 


cm/sec. 
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4. Referring to Figure 3.1, for the case in study where 
2.1 < D < 3, the minimum velocity is approximately 13 
cm/sec and the maximum is 20 cm/sec. 

5. Using the Nomogram (Figure 3.2) for Determining Normal 

Aortic Diameter (Aortic Arch) from P.D. Dr. E. Strehler@° 
the diameter is found to vary from 2.1 cm up to a maximum 
of 3 cm. 

6. Human pulse rate varies from 50 to 100 beats per anes 


Substitution of this data into equations (2.3) and (2.4) 


yields the following ranges of study: 
920 he < 1823 (32) 
9-47 <ta < 193)4 (3:2) 


The flow from the heart will be approximated as a sinusoidal 


type flow as illustrated in Figure 3.3. 
The normalized velocity fluctuation is defined by the 


following expression: 
24 ioe 


where U is mean velocity and U' is the fluctuating component of 


velocity. Kuchar and Scala” report the physiological range for 


dX to be as follows: 
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3.1: Schematic Graph Showing: Broken Line, the Changes 
in Relative Total Cross-Sectional Area (on a 
Logarithmic Scale) of the Vascular Bed; Solid Line, 
the Mean Velocity in the Different Categories of 


Vessel, by A.C. Burton’? 
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Tec he? (3.4) 


Weiting’ found that a 36.7% glycerol-aqueous solution 
was a good hydraulic analogue for blood. The specific gravity of 
this mixture is 1.10 and the viscosity is 4.5 centipoise. Calcu- 
lations of these quantities are in Appendix A. This viscosity falls 
within the range accepted for human blood (3.5 to 5.4 cp.) and the 
specific gravity is reasonably close (1.048 to 1.066). . 
Substituting the kinematic viscosity of the blood analogue 
(v = 0.045/1.1) and the main branch diameter (D = 2.8 cm. = 


1.113") of the experimental model into the relationships for Re and 


a yields the following relationships: 


Re = 185.2 x Q (355) 


(LAS HA (350) 


R 
" 


Where w is measured in cycles per minute and Q is measured in 
liters/min. 
Substituting equations (3.5) and (3.6) into (3.1) and (3.2) 
yields the ranges of flow rate and pulse rate the experimental model 


must be capable of: 


2205 8< 0a 9509 (327,) 


‘demecarinyccoumt Wj he (3.8) 
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In the experiments the required range of the normalized 
velocity fluctuation, A, was covered by the use of different sized 
cams. In Appendix B the interdependence of the length of the piston 
stroke, 1, the velocity fluctuation, X, the flow rate, Q, and the 


pulse rate, w, is outlined. The relationship is: 


1 = 0.770 20/w (3.9) 


where 1 is measured in inches. 


3.2 Description of the Apparatus 


3.2.1 Overall Description of the Apparatus 
Figures 3.4 and 3.5 show the layout and a photograph of 


the apparatus, respectively. The apparatus can be divided into three 


sections: 


1. The flow unit, 
2. The flow rate measurement unit, 


3. The test section. 


Sections (2) and (3) were placed on a table. The test 
section was connected to the flow unit by copper pipes and flexible 
tubings of 1.25 inch internal diameter. 

The inlet and outlet connections of the bifurcation were 


made of surgical tubing.* The upstream branch was approximately 


a ec er streteneesmnennmnntnnesnns nnn teeta eee 


* "Jayon" Surgical Tubing Johnston Industrial Plastics, Edmonton, 
Alberta. 
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twelve feet in length and as straight and horizontal as possible. 


The length ratio of the apparatus was 


Pee) be -90).0038 
as compared to the value of 0.0193 calculated by Dr. Kuchar’. The 
inlet tubing was suspended by twine from a separate stand to appro- 
ximate the physiological case in that a certain degree of Shey 
movement is possible. In addition, the bifurcation was set on 
three rolls and levelled. 

The constant head tanks were made of plexi-glass and manu- 
factured in the machine shop of the Mechanical Engineering Department 
of the U. of A. The tanks were designed and located in such a way 
that the flow lines from the bifurcation to the constant head tanks 
were all approximately of equal length. The weirs of the constant 
head tanks were adjusted to the same constant level (See Figure 3.6) 
to obtain the desired pressure at the bifurcation. It was maintained 
by measuring with a level set.* Pressure upstream of the bifurcation 


inlet was measured using a U-tube manometer. 


* Level set, from the Civil Engineering Department of the U. of A. 
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3.e.c che Flowsunge 
A centrifugal pumps ! located close to the reservoir 

(See Figure 3.4 and 3.5) was connected in series with a rotameter*¢, 
a damping tank (See Figure 3.7) and a pulsation unit. The Calibra- 
tion of the rotameter is shown in Appendix D. The damping tank 
was made of three inch internal diameter plexi-glass tube and a 
height of ten inches. One third of this damper was filled with mix- 
ture, so that the remaining air cushion absorbs the waves from 
the pulsator. A cam-driven piston-cylinder arrangement provided 
the sinusoidal flow. The cylinder was made of four inch I.D. stainless 
steel pipe, both ends of which were threaded to facilitate disassembly 
should any repair or alteration be necessary. The piston was equipped 
with a leather cup seal and the cam follower was a roller-type. To cover 
the range of study, as defined by the dimensionless parameters Re, a, 


and X, various cams were used. 


A variable speed D.C. motor? was used to power the cam. 
To produce the desired range of pulse rate, a 30:1 worm gear reducer*" 


was installed to transmit power from the motor to the cam. Figure 3.8 


x! Centrifugal pump "Jabsco" bronze self-priming motor pump unit, 
model 11810-0003, 1/4 H.P. D.C. motor, Electrical Industries 
Ltd., Edmonton. 

*2 Rotameter, Model 10A 3565A, Max. flow up to 2.7 USGPM, Fischer 

and Porter Ltd., Edmonton, Alberta. 


43 2 H.P., D.C. Motor, Model G.P.-100, General.Electric Ltd., 
Edmonton, Alberta. 


Crofts Model 22 Worm Gear Reducer, Renold Chains Ltd., Edmonton, 
Alberta. 
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is a photograph of the flow unit and Figure 3.9 is close-up view 
of the same unit with a typical cam in place. 

The flow rate, and consequently Re was controlled by 
three valves, two (A and B) located at the outlet of the pump, 


and a third (C) located on the by-pass. 


3.2.3 The Flow Rate Measurement Unit 
The flow was led from the constant head tanks to graduated 
vessels by surgical tubing as shown in Figure 3.10. Flow rates were 
determined in each branch of the aorta by simultaneously closing the 
drains of the measurement tanks and measuring the elapsed time to 


collect a given quantity of fluid (varying from 2 to 6 liters). 


3.2.4 Test Section 

The epoxy resin bifurcation, shown in Figure 3.15 part (H), 
was used for flow visualization. A detailed description of the 
development of the artificial aorta follows in the section 3.4. 

The inlet to the test section was a 1.25 inch I.D. plastic 
pipe while the I.D. of the bifurcation was 1.102 inches (2.8 cm.). 
Each bifurcation branch discharged into a constant head tank (See 
Figure 3.11). The five tanks were designed to approximate the dis- 
charge of fluid from the bifurcation into five infinitely large 
reservoirs. Four reservoirs were of equal size. The fifth tank 
was larger because it was connected to the main branch of the aorta. 
The inlet to each tank was located far from the sidewalls in order to 


avoid wall effects. 
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3.3 The Experimental Fluid 

Water was chosen as the experimental fluid for the first 
part of the experimental work where qualitative observations were 
made of the flow characteristics in the bifurcation. 

For decreasing the surface tension an agent*! was added 
to the water. This was done to ease the flow over a weir, between 
internal walls of the bifurcation and constant head tanks made of 
plexi-glass. The viscosity was measured with a Cannon-Fenske- 
viscometer*® atezZene 

For the quantitative experimental studies a mixture of 
water and glycerol was chosen as the experimental fluid. To simulate 
more closely real blood flow, the author adopted Weiting's?* recom- 


mendation that a 36.7% glycerol-aqueous solution was a good hydraulic 


analogue for blood. 


3.4 Description of the Development of the Artificial Aorta 


A sketch of the aortic arch, giving the relevant dimensions, 
was obtained from Dr. N.R. Kuchar. Further assistance was received 
from Dr. Kocandrle, Dr. Kubac and Dr. R.D. Laurenson. As a result 
of this assistance an artificial aortic arch was designed (See 
Figure 3.12). This design created a basis for the preparation of 
a model. Figure 3.13 shows the basic steps in the construction of the 
aortic arch. 

Part A of Figure 3.13 shows the basic part machined to the 
dimensions of inlet and main outlet branch of the aortic arch. The part B, 
*1 Photo-Flo 200 Kodak, Edmonton Photo Supply, Edmonton, Alberta 


*2 Cannon-Fenske viscometer No. 50J 657, Chemical and Petroleum 
Engineering of the U. of A. 
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was machined out of aluminum. 

The two halves of the model shown in part C were obtained 
using part B. Figure 3.14 shows part C in detail. The spaces of 
the cutting were filled with soldering tin. 

The box shown in part D was designed to enable easy re- 
moval of the epoxy-resin cast aorta shown in part E. The cast aorta, 
as shown in part E, was used for production of several paraffin models, 
part F. These models were repeatedly dipped into the Apex Latex* and 
dried so that subsequent layers built up a flexible wall of an aortic 
arch. After drying the paraffin from inside was melted out in hot 
water. This way a final product of flexible aorta (part G) was 
obtained with the same inside dimensions as part E. 

Part E was worked further by milling, grinding and polishing. 
The two half bifurcations were finally bolted together and/or sealed 
with chloroform to form part H, the final test section. This was 
transparent and used in the visualization experiments. All final 
types are shown in Figure 3.15. The glass model**(Part I) was used 
in preliminary studies. 

Unfortunately a flexible and transparent model could not 
be developed because of lack of suitable materials. Roneved: it 
would be desirable in further studies in this field to develop such 
an aortic arch. 

Generally, the final test model (H) is slightly different 
from the actual human aorta because all branches are in the same plane, 


whereas in the human body they are declining in various small angles. 


rear er eo eee 
* Apex Latex, supplied by Continental Rubber Co. of Canada Ltd., Edmonton, Alberta 
**Glass model was made by Glass Shops at U. of A. 
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3.5 Flow Visualization 
Resin particles and the "Hydrogen bubble technique" were used by 
previous investigators for flow visualization (Rodkiewicz - Howell and 
Rodkiewicz - Roussel). The hydrogen bubble technique appears to be much 
more efficient. The disadvantage of this method is that water must be used 
as the experimental fluid. If a water-glycerine mixture is used, the con- 
ductance, and therefore, the production of bubbles decreases. 
A tungsten wire of 0.002 inch diameter was used as a cathode and 
a copper ring was used as an anode. The following equations describe the 
reactions which take place. The oxidation potentials of tungsten and copper 
are +.12V and -.34V, respectively, 
Anode: Cu + 20H > Cu (OH). +2e 
Cathode: 2H" + 2e +H, (gas) 


Water equilibrium: H,0 # HeeOne 


The hydrogen bubble technique utilizes the electrolysis of water 
to introduce hydrogen into the flow. These bubbles, formed at the fine-wire 
cathode are swept off by the flow. The bubbles produced were 0.001 inches in 
diameter and have a maximum vertical velocity, due to buoyancy effects, of 
0.0011 fps and reach a velocity of 98 percent of the fluid velocity within 1 ms, 


25 according to Rodkiewicz - Roussel. 


as explained by W. Davis and R.W. Fox 
The wire has, due to its size, an insignificantinfluence on the flow 
itself. This constitutes one of the best advantages of this method, which also 
provides a great flexibility to yield qualitative and quantitative results. 
Qualitative results can be obtained by continuously producing bubbles. 
The success of this visualization method is highly dependent on 


correct illumination. A sharply collimated sheet of powerful light and dark 


backgrounds are necessary. Other light sources must be suppressed. 
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The tungsten wire and the copper ring were built-in to a special 


electrode housing (Figure 3.16). The electrode housing was designed so that: 
1. The wire could be replaced. 
2. It could rotate about its axis in order to generate 
sheets of bubbles in any of the diametrical planes. 
The electrode housing was located as shown in Figure 3.11. 
The generation unit, which allows the operator to have 
complete and precise control of the frequency and the duration of 
the pulse producing the bubbles, is shown in Figures 3.17 and 3.18. 
From this generation unit the electrolysis DC current is brought to 
these electrodes. The main part of this generation unit is the 
constant DC potential power supply*!, 
When the switch is turned on by the pulse generator*@ which 


3 360 volts are applied be- 


is driven by the Wide Range Oscillator* 
tween the two electrodes. The oscillator regulates the frequency of 
the pulses, while the pulse unit controls their duration. A reversing 
Switch was installed so that the wire could be cleaned. 

The photographs of the steady flow streamlines were taken 


with a reflex camera*" loaded with high speed film (Kodak Tri-X Pan, 


ASA 400). A movie camera*? (Kodak 16 mm, Tri-X, Reversal film, Type 7278) 


x! Lambda Regulated Power Supply, Model C 882M, Lambda Electronics Corp. 
xe Pulse Generator, Type No. 1217 C, General Radio Company. 
*3 Wide Range Oscillator, Model 200 CDR, Hewlett Packard Ltd. 


«4 Canon Camera, Model F Tb, f 1.4; 50 mm. 


*° Movie Camera, Type: Bolex, H 16 Reflex, 10 mm - wide angle Bolex 
Lens. 
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was used for the pulsatile flow studies. 
One slide projector (total power 450 W) and one car head- 
light were used to give a collimated sheet of light in the same 


plane as the bubble sheet. 
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CHAPTER IV 
EXPERIMENTS, RESULTS AND DISCUSSION 


4.1 Experimental Procedure 
4.1.1 Preliminary Procedure 


A desirable inlet pressure in the bifurcation was maintained 
by the constant head tanks. The liquid column height was measured with 
the "level set" instrument and maintained at the equivalent of 115 mm Hg 
for all measurements. 

At the beginning of each experimental session it was neces- 
sary to go through the following procedure: 

1. The circulation pump was turned on to circulate the 
mixture of glycerine and water for at least 30 minutes 
in order to obtain a homogeneous medium. 

2. The dumping tank was filled to approximately 1/3 full 
with the mixture (2/3 left for air). 

3. The density of the blood analogue was checked with the 


hydrometer and adjusted as required. 


4.1.2 Actual Experimental Procedure 


1. The appropriate cam was installed. 
2. The centrifugal pump was started and the required 
flow rate was set on the rotameter using valves A 
and B (See Figure 3.4). Valve B was set approximately 


and the needle valve A was used for fine adjustment. 
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3. The pulsatile unit was set at the desired rotational 
Speed and checked with a stop watch. The de-aeration 
of the pulsator was done using the valve on the top of 
the pulsator. 

4. The drains from the calibrated tanks were closed 
Simultaneously and the elapsed times to gather two 
kiters in tanks 1, 2, 3, 4 andsix liters for tank 5 
were measured by using five stop watches. 

5. The temperature of the solution was measured in the 
main reservoir. The temperature remained constant at 
76°F for the duration of the measurements. 

6. The flow rates in each branch Qy> Qo» Q3> Qy> Qc» the 
mass flow ratio Yy> Yo Y3> Yq> Y5 and total flow Q 
were computed. In general there was at most a + 2% 
difference between the measured flow rate and the flow 


rate given by the rotameter. 


4.2 Photographic Procedure 


As mentioned in section 3.2.2 water was used as the expe- 


rimental fluid to obtain the photoaraphs of the distribution of the 
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flow through the artificial ascending aortic arch. The conditions 
under the photographs were taken are as follows: 
1. Only one type of cam, of length, 1 = 0.27 inches 
was used. 
2. Before taking the photographs, the Re was determined 
(See section 4.1.2 of this chapter). 
3. All photographs were taken with the "Cannon" Reflex 


Camera of the "Bolex" Movie Camera. 


4.3 Experimental Results 
All of the experimental data are listed in Appendix C. 


4.3.1 Steady Flow 
All the data has been taken using a mixture of glycerine and 
water however, photographs were obtained with water alone. The 
bifurcation inlet pressure was 115 mm Hg. The experimental values were 


used in plotting Figure 4.1 listed in Table C-4a. 
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4.3.2 Mass Flow Ratios as a Function of the Unsteadiness 
Parameter and the Reynolds Number 


The normalized fluctuation component of velocity, A, was 
kept constant at a value of one. 

The dependence of ¥1> Yoo Y3> Yq> and Ys On a was examined 
at 


Re = 800, 1300, and 1800. 


The size of cams used (different length of stroke, 1), w, a, 
and the experimental values obtained for 1) You ec atae and Ys are 
given in Tables C-4.2, C-4.3, and C-4.4. The data of the reproducibility 
tests are marked with a comma as illustrated in Figure 4.4. 

Figure 4.3 shows the relation between Y¥y> Yoo Y3> Yq» and 
Me and a. The Reynolds Number was held constant at 1300 and A was 


held constant at 1.5. The experimental data are given in Table C-4.1] 


4.3.3. Mass Flow Ratios as a Function of the Reynolds Number 
and Unsteadiness Parameter 


The normalized fluctuation component of velocity, A was kept 
constant at a value of two. 

The dependence of Yy> Yoo Y3> Yae and Ys on the Reynolds 
Number was examined at a = 10, 15, and 19. The data in Tables C-4.7, 


C-4.9 was plotted in Figures 4.4, 4.5 and 4.6. 


Figure 4.4 shows the variation of vy (mass flow ratio 


associated with right subclavian) with Re as a function of a. 
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Figure 4.2: Dependency of Mass Flow Ratios on a; 
Rev= (300,07 = 15; 


15) 
15 


8 10 12 14 16 18 20 
O Re = 800 o Re = 1300 A. Re = 1800 


Figure 4.3: Dependency of Mass Flow Ratios on a3 
Re = 800, 1300 and 1800, A = 1. 
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Figure 4.5 shows the variation of Yo (mass flow ratio 
associated with the right carotid) with Re as a function of a. 

Figure 4.6 represents the experimental results for Re and 
Y3 (mass flow ratio associated with the left carotid). 

Figure 4.7 shows the variation of Yqs (mass flow ratio 
associated with the left subclavian) with Re as a function of a. 

Figure 3.8 gives the variation of Y5» (main mass flow ratio 


associated with the thoracic aorta) with Re, as a function of a. 


4.3.4 The Mass Flow Ratios as a Function of the Reynolds 
Number, and the Normalized Fluctuating Component of 
Velocity 


The aera dances parameter was kept constant at a = 10. 
A detailed description of the mass flow ratio and their relation to 
the branches of the bifurcation is found in section 4.3.3. The 
dependence of ¥1° Yoo Y3° Ya> and Ys on Re was examined at A = 1, 
1.5 and 2. The experimental data are compiled in Tables C-4.5, 
C-4.6 and C-4.7. 

Figures 4.9, 4.10, 4.11, 4.12 and 4.13 show the relation 


of Y1> Yoo ¥3> Yq and ‘its to the Re and dX for a = 10. 


4.3.5 Description of the Flow 


Similar characteristics of the flow in an artificial aortic 


arch can be found in steady flow as well as in the forward flow of 
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the sinusoidal pulsatile flow, as illustrated in Figure 4.14. 

Four separation regions are shown in this sketch. One 
is in the area where the right common carotid artery divides. Another 
clearly described region is seen in the area of the innominate artery 
at the junction of the right subclavian (Branch 1) and the right 
carotid. Thin separations also occurs in the left carotid and the left 
subclavian (Branch 3 and 4). These regions, however, though difficult 
to sight, could not be photographed due to light reflection at the 
branch edges. 

The flow pattern of the back pulsatile flow is shown 
in Figure 4.15. In this situation the separations regions seen in 
the steady and forward pulsatile flow, do not appear. However, a 
separation region does occur On the internal edge close to right 


innominate artery. 


The dependence of the steady flow characteristics on 

_ different Reynolds Number, is shown in Figures 4.16, 4.17, 4.18, and 
4.19. From Figure 4.20 one can infer the distribution of the 
streamlines in the aortic arch by the steady flow. The flow 
distribution in one complete pulse of the pulsatile flow is shown 


mn) Ejgure 421. 
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Figure 4.14 Flow Pattern for Steady and Forward Flow by the 
Pulsatile Flow. 
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The movie pictures in Figures 4.20 and 4.21 were taken 
at approximately the same Reynolds Numbers, (Re = 1230, Re = 1180) 
and show distinctively the different behaviour of steady and pulsa- 
tile flow. 

The appearance of the double-helicoidal flow in the main 
branch down stream of the lower aortic arch at different Reynolds 


Numbers is shown in Figure 4.22. 


4.4 Observations and Discussions 
4.4.1 Mass Flow Ratio 
The study of the diagrams given in Figures 4.1 through 4.13 
leads to the following observations: 


a. Mass Flow Ratio, vy 


1. Figures 4.2 and 4.3 indicate that Y) is independent 
of a for. a given Re, in range of 10 < a < 20, and 
decreasing function of a and Re in lowest range of a 

(oe< One C)einaG 1S; YY slightly increases as a 
decreases and vy moderately increases as Re decreases, 

in this range. 

2. In the particular case of sinusoidal flow through 

the aortic arch vy is a decreasing function of Re (Figure 
4.5), fora given a. That signifies if a is increased 


and X is kept constant. 
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Figure 47.12% 


Flow Pattern for Back Flow by the Pulsatile Flow. 
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Fiqure 4.16 Steady Flow Characteristic or \ 
Reynolds Number 


c) Re = 1570 d) Re = 1670 


Figure 4.17 Steady Flow Characteristics for Various 
Reynolds Number 
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Figure 4.18 Steady Flow Characteristics for Various 
Reynolds Number 
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Fiqure 4.19 Steady Flow Characteris ’ 
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3. It was found (Figure 4.10) that ‘A is a decreasing 
function of Re and the dependence of Y; on A is such 
that vy decreases with increastng A, for a = constant, 


up to Re = 1500. For higher Re, Y) is a decreasing 


function of Re only. 


4 For steady flow ‘a is a slowly decreasing function 


of Re, (Figure 4.1). 


b. Mass Flow Ratio, Y9 


The dependence of Yo» are Similar to those of yy except 
as follows: 

1. For a sinusoidal flow the dependence Yo on Re 

(Figure 4.5) is much stronger than the dependence of 

Yon Re (Figure 4.4). This may be caused by the fact 

that the separation region created at the junction of the 

innominate is forcing the flow to the right carotid, 

(Branch 2). | 

2. Figure 4.10 indicates that Yo is a decreasing function 

of Re. Above a certain Re, yo is a decreasing 

function with increasing A. The dependence of Yo on r 


is not as strong as that of ¥, on dX as shown in Figure 


4.9. 


c. Mass Flow Ratio, ¥3 


i 


The same dependence can apply as in section (a). 
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d. Mass Flow Ratio, Yq 


—_—_——— 


The same situation occurs as in section (a) except for the 
diagram, Yq versus Re, (Figure 4.12). Here Yq is not so clearly a dec- 
reasing function of A as was the case for Y]» Yg2 and y,. The dependence of 


Yq on A is not apparent and cannot be determined without further testing. 


e. Mass Flow Ratio, Ys 


1. Figures 4.2 and 4.3 indicate that Ys is not a function 

of a for a given Re. 

2. For the sinusoidal flow through the main branch of 

the artificial aorta, fs is not a function of a, and 

v5 is a function of Re only. The velocity fluctuation 

being kept constant. 

3. Figure 4.13 indicates that Ve is not a decreasing 

function of as was the case for Yy> Yoo and Y3- The 

dependence of Ys on X is not apparent and cannot be 
determined without further testing. 

4. In steady flow, 5 is a function of Re. 

In previous sections the results were described and the 
interpretation of the physical meaning can be drawn in the following 
way as indicated by Dr. Kuchar*. 

When the geometry of the bifurcation is a highly curved 
one, the effect of Dean number, K, will increase with increasing Re. 


I 


* Private communication between Dr. Rodkiewicz and Dr. Kuchar. 
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As a result, the centrifugal force will aid in driving blood into 
the head and arms. The pressure gradient at the centre zone of the 
aortic arch will create secondary flow from the internal to the 
external wall of the aortic arch. This occurs where the first three 
branches (innominate aorta, left carotid, and left subclavian)are 
located. 

This would indicate that the flow ratios of the first four 
branches (the upper extremities) decrease with increasing Re, at the 
expense of the flow ratio of the main branch. In summary it appears 
that as the Reynolds Number increases the mass flow ratios of flow to 
the first four branches decrease and the flow ratio of the main branch 
increases correspondingly. That is, as the total blood flow is 


increased a smaller portion is diverted to the upper extremities. 


4.4.2 Separation Regions 
a. Separation Created by the Steady Flow and Forward 


Pulsatile Flow 
In this particular case, the photographs and observations 
clearly prove the existence of four interdependent separation regions: 
1. The largest separation region is in area of internal 
wall of upper and lower aortic arch and continues into 
or ends in the vicinity of the upper thoracic aorta. 


The termination of this region is difficult to determine. 
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It can be seen from Figures 4.17 through 4.21 that 

by increasing the Reynolds Number, point X of the 
separation (Figure 4.15) is moved further along the 
internal wall upstream and its thickness Te becomes 
greater. The cross-sectional shape of the separation 
region could not be determined and the shape shown has 
been assumed. Observation of the secondary flow how- 
ever seems to indicate a profile of this shape. 

2. The other recognizable separation occurs in the area 
of junction of the innominate. From the photographs men- 
tioned above it can be found that by increasing the 
Reynolds Number the point Y of the separation moves 
further along the wall of the innominate upstream, and 
the thickness T19 becomes larger. The tail of the sepa- 
ration moves the same as point Y, to the upstream with 


increasing Re. 


3. The two smaller separation regions appear in branches three 
and four on the right walls of these branches (Figure 4.14 and 

4.15). The separation thicknesses T3 and Tq are so small] that 

no measurement was possible. In general only qualitative 


results could be obtained due to smallness of the regions. 


b. Separation Created by Back Flow of the Pulsatile 


Flow 
During back flow all the previous described separations 
disappear and a new separation region is created as Shown in Figures 


4.16 and 4.22. It was found that by increasing Reynolds Number the 
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Point of separation Z is moved downstream and the separation region 
became more convex. 

A vortex formation was seen in the area at the end of the 
main branch. Its occurrence can be seen tn Figure 4.21. 

Inside the main branch and branches of the left carotid 
and left subclavian, the double-helicoidal flow was observed. In 
the area downstream of the lower aortic arch, as the flow entered the head 
tank, photographs were taken perpendicular to the cross-section of 
the main branch. In this particular cross-section the axis of double- 
helicoidal flow is oriented as shown in Figure 4.22 or as sketched 
in Figure 4.14 (Section E-E). This kind of flow is typical of flows 


in curved pipes as explained by H. Schlichting !®, 


4.5 Limitation of the Apparatus 


Several test runs were repeated to check if the test runs 
were reproducible. The results were as follows: 

The results of the check of the experimental sensitivity 
were found to agree to within Eee he 

The reproducibility is directly proportional to the care 
taken to obtain precisely the same level for all five constant head 

tanks, because the flow through the bifurcation is highly dependent 

| on the resistance in each branch. 

The visualization technique becomes useless below Re = 600, 


due to the growing importance of the buoyancy effect on the hydrogen 


bubbles. 
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CHAPTER Y 
CONCLUSIONS 


5.1 Summary 


For the physiological range of unsteadiness parameter, 


Reynolds Number, and for normalized velocity fluctuation, the experi- 


mental findings can be summarized as follows: 


LE 


The mass flow ratio, yy] iss 

a. Not a function of a for a given Re in range of 
TOg<fam<7 20: 

b. Inversely proportional to a in the lowest range - 
(8 < a < 10)for a given Re. 

c. A decreasing function of both Re and i, for 

a = constant up to the range of Re = 1500. For 
increased Re, Y] is a decreasing function of Re only. 
d. In steady flow, a decreasing function of Re. 

The mass flow ratio, Yo has same dependences as yy 
except as follows: 

a. The dependency of Yo on Re iS more pronounced. 
This might be caused by the separation created in 
region of junction of the innominate. 

b. The dependency on id, is not so strongly indicated. 


The mass flow ratio, Y3 has the same dependences as in 


case of Yq: 
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The mass flow ratio Yq has the same dependences as 

Y) except the dependence on id is not apparent and can- 
not be determined without further testing. 

Results of the other mass flow ratios indicate a 
dependence on A, however Figure 4.13 shows no discern- 
able trend; perhaps due to the limitation of the expe- 
rimental sensitivity. It is shown, however, that Ys 

is a function of Reynolds Number. The mass flow ratio 
Ys 1S not a function of a at least in range 10 < a < 20 
for a given Re. 

There exist four interdependent separation regions 

in the steady and forward pulsatile flows. 

With back pulsatile flow, all separations mentioned above 
seem to disappear and are replaced by one located at the 


outside wall just before the innominate branch. 


A visualization technique proved the vortex formation 
observed in the separated region of the main branch 
under lower aortic arch. 

The double helicoidal flow has been observed in branches 


4 and 5 in steady and forward pulsatile flows. 


5.2 Suggestions for Further Work 


1h 


It would be of interest to obtain clearer photographs 
and movie pictures so that better visualization of the 
distribution of the flow could be observed. It might 


be interesting to magnify proportionally the artificial 
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aorta and improve the illumination (for instance from 
all four sides) to remove reflected light from the side 
walls branches. 

Repeat the same tests with the artificial flexible aorta 
(a flexible aorta was developed by the author with the 
assistance of Mr. Golls) and results compared with those 
obtained in this thesis. 

Repeat the tests performed using redesigned aortas with 
different geometry, mainly with different curvature of 
aortic arch to study how the curvature would influence 
the dependence of mass flow to the branches, (Dean number, 
K; or Modified Dean number, Km, versus mass flow). 

Tests should be performed which will better simulate 
conditions for "real" flow through the aorta. This data 
obtained can then be compared with results obtained from 
sinusoidal flow. 

Finally simulate the experiments "in vitro" 

and compare the results with the artificial one. To 
simulate the "real" flow in aorta, the author suggests 
using the same constant head tanks, graduated vessels 
and the reservoir, as used in this experiment. The 
artificial left ventricle with its artificial valves 

put in series as follows: constant head tanks from 
overflow to graduate vessels - reservoir - the free 

flow lead into the artificial left ventricle from aorta 


back into the bottom of constant head tanks. 
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APPENDIX A 


APPENDIX A 
SPECIFIC GRAVITY AND VISCOSITY OF 
THE BLOOD ANALOGUE 


This appendix is a condensation of the detailed studies 


done by Rodkiewicz-Howell !¢, from which the following physical 


constants have been obtained: 


Tk 


The blood analogue used was a mixture of glycerine 
and water of which 36.7% was glycerine by volume. 
The specific gravity of glycerine, measured with 
an Exax Heavy Liquid Hydrometer, was 1.270. 


The specific gravity of the mixture is as follows: 


.367 (1.270) + .633 (1.000) 


466 + .633 


1099 


The weight percent of glycerine of mixture 


Thus 42.4 1b. of the glycerine and 57.6 1b. of the water 


86 


were weighted for each 100 1b. mixture. The value 1.1 of 


the mixture was checked with the Exax Heavy Liquid 


Hydrometer before each experiment was carried out. This 


mixture gives a viscosity of the blood analogue of 4.5 


centipoise. 
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APPENDIX B 


PISTON STROKE CALCULATIONS 


The length of the piston stroke was varied by using cams 
of different sizes, so that the effect of changes in the normalized 
velocity fluctuation could be studied. The fluctuation component, 
A is also dependent on the pulse rate, w and the flow rate, Q. 

The value of A will change if w or Q is changed, even 
though the cam is not changed. 

The relationship between Q, w, A, and the length of 
stroke, 1 is determined as follows: 

Consider the instantaneous flow rate simulated on Figure 


3.3. It can be described as equal to 


Were as 
Q + a~ Sin (7 t) (B.1) 
0) 
where: Q' is fluctuating component of the flow rate, 
ee is the time required to complete one pulse, 
t is the time variable. 


Since Q and Q' are directly proportional to the velocity, 


U, and fluctuating component of velocity, U', it is possible to write 


so that 
Or eAQ. (B.2) 
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Therefore, the instantaneous flow rate becomes 
Q [1 + x sin (=) al 
fe) 
The required swept volume of the piston equals the shaded 


area Of Figure 3.3. 


The required swept volume becomes 


Tg 
2 
XQ sin (25) t dt 
Z if 
O 
To 
2 


£1 Be ED + 32 72 cos 

eet 0 0, ag 70 

bape 

“nat, (B.3) 
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The time required for one pulse equals the reciprocal of 
the pulse rate. 


Therefore 


=] 
Tyts (B.4) 
Substituting equation (B.4) into equation (B.3), the 
required swept volume becomes 
ite (B.5) 


The required length of stroke is obtained when the volume 
is divided by the cross-sectional area of the piston (diameter = 4 
inches and applying the appropriate conversion factors. 


That is, 


oe x 61.025 (inch)? 
ie 
T (2)° (inch)° 


Bao 
1 = 0.77 “2 (inch) (B.6) 
The minimum 1 was found to be 
x ] 2.85 
Ves VEN Seng: 
1_.. = 0.019 inches, (B.7) 
min 


(2,8) 
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and the maximum 1 was found to be 


] ee 0.541 inches. 
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(B.8) 
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Table C-4.a 


Dependence of Mass Flow Ratios on Re for Steady Flow 
of Mixture (Glycerine + Water) at a Static Pressure of 115 mm Hg 


Ss yy 12 ¥3 Yq Y5 
370.4 129 196 157 196 333 
463.0 132 181 155 190 356 
655.6 127 169 145 184 373 
648.2 123 158 142 173 389 
740.8 126 AS 138 176 410 
833.4 122 150 131 170 419 
926.0 120 147 130 168 425 

1018.6 120 143 128 ce 439 
1111.2 17 141 122 168 448 
1203.8 123 eo 129 163 456 
1296.4 7 135 126 164 463 
1389.0 114 133 125 161 476 
1481.6 121 126 119 159 483 
1574.2 114 130 aie 56 470 
1666.8 116 128 119 157 476 
1759.4 12 124 15 159 485 
1792.7 15 125 114 155 495 


Re = 185.2 x Q; Q (liter/min.) 
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Table C-4.1 


Dependence of Mass Flow Ratios on 
a for Re = 1300 and A = 1.5 


O235 23 8.58 .118 aey, 124 
xOR35 Zara 8.58 ma F/ -136 .126 
OF33 24.51 8.87 114 ise) wi22 
0.31 26.09 9.12 114 i133 SZ 
OR] 29.96 9.78 114 Bi Sy) .119 
*0.27 29.96 9.78 112 133 121 
0225 32.36 10.15 114 “133 ah 
Oe73 S517 10.58 111 130 116 
*0.23 a5 e7 10.58 112 130 120 
0.20 40.45 1137 113 128 122 
0.18 44.94 11.97 112 131 118 
0.15 53.93 3e10 i 130 122 
*0.15 53.93 13.10 112 129 120 
O22 67 .41 14.67 114 133 122 
0.10 80.90 16.06 111 135 123 
0.08 101.12 17493 114 132 123 
x0 .08 LOZ 17.93 114 13] 122 
0.07 115.57 19.15 113 133 27 


Q = 7.01 liter/min. 
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Table C-4.2 


Dependence of Mass Flow Ratios on 


We 


Q = 4.31 liter/min. 


a for Re 


800 and A = 1 
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Table C-4.3 


Dependence of Mass Flow Ratios on 


a for Re = 1300 and A = 1 


£03 


Q = 7.01 liter/min. 
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Table C-4.4 


Dependence of Mass Flow Ratios on 


a for Re = 1800 and A = 1 


Leh: 
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Q = 9.71 liter/min. 


Table C-4.5 
Dependence of Mass Flow Ratios on Re 


for a = 10 and X= 1 


l(ins) Q(1.p.m.) = Re Yy V9 Y3 
O23 EY! 1736 107 .126 oERZ 
*0 23 9337 1736 109 Shen ~LT6 
0.20 8715 1509.5 TZ | Si ie 
0.18 dian 1358 Byte 13] -\2] 
vi airs 7*833 1358 Pa a) pics ace 
O15 6.11 I32 Pa AL: .136 ee 
0.14 5270 1056 io 14) .124 
Onts 5.39 981 116 . 140 a25 
Ost 4.88 905 si llee . 144 oe 
OF Wl 4.48 830 .118 .149 eW29 
Ore td 4.48 830 ai! . 148 loo 
0.10 4.07 ie .124 Al sie: 133 
=O 10 4.07 15 119 .150 ioe 
0.09 3.66 679 .128 . 159 139 
=0 209 3.66 679 WAS) wiv) Bey! 
0.08 Saco 604 2 eh . 160 lec 
“O05 Seed 604 aliZs . 164 .140 
0.07 2450 528 .129 . 169 .144 
$0. 07 apa 528 us ae ~145 


w = 31.38 revolutions/min. 
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Table C-4.6 


Dependence of Mass Flow Ratios on Re 


for a = 10 and X = 1.5 


Geo0 TO%32 1912 
0.36 9.78 1811.4 
0235 9.50 176] 
0.34 9.24 Toe? 
0233 8.96 1660 
0.32 8.69 Toro 
O33) 8.42 1560 
O32? Tva3 1358.5 
GeZ25 6.79 1258 
*0225 6.79 1258 
0823 6.24 aber 
*0k23 6.24 T157 
0.20 Ges 1006 
0.18 4.89 906 
*0.18 4.89 906 
Jag 4.07 inate 
*0.15 4.07 755 
0.14 3.80 704 
*0.14 3.80 704 
Ded je 3705 654 
*0sl3 3753 654 
O2tZ Shela: 604 
0.11 2.98 Sod 
*0.11 2.98 aya. 
0.10 (Dros 503 
*0.10 Ae 7ih 503 
w = 31.38 revolutions/min. 
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Table C-4.7 


a= 10 and A = 2 
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Tins) .0(1.p.m) Re yy Yo Y3 v4 Ys 
0.47 OF 57 1773 .6 .108 AWS e120 GS .482 
0.46 Or37 1736 «143 .129 ss . 166 .467 
0.45 Oe 1698.2 .110 eico silts .168 477 
0.44 8.96 1660 Bil -123 el Roe -470 
0.43 8.76 1622.7 5110 .128 .118 .164 -472 
0.42 oe 1585 aa A esl mics . 164 .474 
0.40 SaL5 1509 ells Pio .120 163 -465 
0.38 7.74 1434 = ae 715" 2121 . 166 -464 
0.35 Fos 1320 16 2138 -120 PA -467 
0.33 Gaiic 1245 P.O witehs) Pliee LOy, -469 
0.32 B02 1207 .6 2 158 ehe2 VAG .464 
0.31 Oro 1170 ou ais) ‘133 .119 . 166 .458 
02% 5.00 1019 ral 2186 lel 67 -436 
0225 5.09 943 ite ehZo ee ie .436 

vad A 5.09 943 114 SUSIE aLe6 . 166 44] 
O23 4.68 868 . 107 a1o7 S122 AG .433 
0.18 3.06 679 Ail .148 -130 175 .409 
Oe15 305 566 .110 By ste: rio .184 .402 

~O. 15 Shes) 566 .114 154 aga .180 .413 
0.14 L100 528 e113 .164 A eye .189 412 

*0.14 2205 528 mbes) sen 27 .188 413 
One 2.64 490 =123 RIA 144 . 183 2304 

xO be 2.64 490 A . 166 .146 mers) .384 


w = 31.38 revolutions/min. 
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Table C-4.8 
Dependence of Mass Flow Ratios on Re 


for a = 15 and A = 2 


1(ins) Q(1.p.m) Re vy Y5 Y3 Hp Y5 
0.20 emell, 1698 .107 pe .114 Eto -484 
mURcO ea 7 1698 alla) ee 6d Br 7 b5S8 49] 
0.18 Sharale) 1528 lal eg 116 .161 478 
20018 Ono 1528 oT . 188 .119 elie ~475 
Ott 6.87 1274 . 109 gl ES YZ 124 i Kot .457 
Ufa) 5.96 1104 iis) a oie) .126 . 166 .470 
Wi ba 520 1019 A fits) as ee) i as Pos .43] 
0.11 5.04 934 val 4) .142 fi32 . 168 -440 
FORK 5.04 934 Aisle 142 nlZo A 434 
0.10 4.58 849 Rg, 142 nIZO . 166 433 
beal (0) 4.58 849 -120 144 129 lee .433 
0.09 4.12 764 Bice 150 .134 174 BN i) 
*0.09 4.12 764 eB aa We) .134 ITS 418 
0.08 3.66 679 124 156 139 176 393 
*0.08 3.66 679 125 159 140 175 390 
0.07 Seal 594 126 163 iy, 180 370 
SUcUy eral 594 .124 . 163 . 140 .182 .376 
0.06 aed fe) 509 5130 sy Ld .148 . 183 .364 
0200 Geto 509 ns bef) oh GES: soe .184 70/0 
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w = 70.62 revolutions/min. 


Table C-4.9 
Dependence of Mass Flow Ratios on Re 


for a = 19 and A = 2 


l(ins) Q(1.p.m) Re vy Y5 73 Y4 
0.14 10.30 1907 110 119 114 15 
Oars 00 1771 114 120 Wb 153 
Or 8.82 1635 Ws 127 116 157 
Oni 8.09 1498 ALLO) 126 117 153 
0.10 730 1361] 114 128 15 160 
xO nL0 Techs) 1361 lO Bilas) lit 161 
0.09 6.62 1226 .118 134 124 . 164 
e209 6.62 1226 Sy ie 134 eli2s . 168 
0.08 5.88 1090 2120 pete: Bee - 164 
*0.08 5.88 1090 e116 Bo al23 . 168 
0.07 Dao | 954 sl20 141 “28 . 180 
0°07 Bato 954 ars log Aa mel Ar 
0.06 4.4] 817 SE 149 e134 .170 
*0 .06 4.4] 817 S141 . 148 cio .170 
0.05 3.6/ 681 117 1 131 178 
=0,05 3207 681 sl22 sou 29 SUAS 
0.04 2.94 545 .126 iz 146 181 
*0.04 2.94 545 126 Iv2 15a 187 
0.03 .20 409 p33 193 ey) els) 
+0703 2.20 409 lied rlles}s) n25.0 9) 


w = 113.3 revolutions/min. 
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APPENDIX D 
ROTAMETER CALIBRATION 


The calibration results for the rotameter are shown in 
Figure D.1. The curve obtained is in linear; the volume rates for 


intermediate rotameter readings are obtained by interpolation. 
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